N95- 13926 


SQUEEZING GENERATED BY A NONLINEAR MASTER 
EQUATION AND BY AMPLIFYING-DISSIPATIVE 

HAMILTONIANS 


V.V. Dodonov 

Moscow Institute of Physics and Technology 
16, Gagarin street, If 0160 Zhukovskiy, Moscow region , Russia 

M.A. Marchiolli, S.S. Mizrahi and M.H.Y. Moussa 
Departamento de Fisica , Universidade Federal de Sao Carlos 
Via Washington Luiz km 235, 13565-905 , Sao Paulo, SP, Brasil 


Abstract 

In the first part of this contribution we show that the master equation derived from the 
generalized version of the nonlinear Doebner-Gloldin equation leads to the squeezing of one 
of the quadratures. In the second part we consider two familiar Hamiltonians, the Bate- 
man - Caldirola - Kanai and the optical parametric oscillator; going back to their classical 
Lagrangian form we introduce a stochastic force and a dissipative factor. From this new 
Lagrangian we obtain a modified Hamiltonian that treats adequately the simultaneous am- 
plification and dissipation phenomena, presenting squeezing, too. 


1 The nonlinear master equation 

Our search for a new class of equations was inspired by the model proposed recently by Doebner 
and Goldin [1], They looked for the most general Schrodinger equation compatible with the 
Fokker- Planck equation for the probability density p(x,f) = | | 2 , 


Pt + V • j = D V 2 p 

(where D is a constant positive diffusion coefficient), and derived the nonlinear equation 


( 1 ) 


ihdxl>(~x.,t)/dt = 
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[~(ft 2 /2m)V 2 + V(x)]V>(x,t) 
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v 2 + | vy>(x,o p 
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rp(x,t) 
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for a particle with mass m moving in a scalar potential V(x). The advantage of eq. (2), com- 
paratively to other ones [2], is its group theoretical origin, the nonlinear term was not simply 
added to the usual Schrodinger equation ad hoc , but its structure was derived from the analysis 
of representations of the Diff[ R 3 ) group, proposed as a "universal quantum kinematical group” 
[1]. The only drawback of eq. (2) is its limited range of applicability, since it can be used only in 
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the coordinate representation. It is desirable to have a more general equation valid in arbitrary 
representations. To obtain such a generalization (heuristically), first, we remove the Cartesian 
coordinates dependence, introducing an arbitrary set of states | z), which form a complete system 
with respect to some measure dp(z), i.e., 


J I z)(z | dfi(z) = 1. 


(3) 


Secondly, we replace the two V - operators in the additional term of eq. (2) by two arbitrary 
(linear) operators A and B. By this way we arrive at the equation, whose nonlinear part, in 
general, is neither Hermitian, nor anti-Hermitian, 


ihj t (z\4>) = (z\H\iP) + tnD 


( z I V») 


(tH £ 1 z H z 1 * I VO 
(V> I z)(z I V>) 


{z\BA\xl>) 


(4) 


Multiplying both sides of this equation by (ip \ z) and integrating over dp(z) with account of (3), 
one can check that the normalization of the wave function is preserved in time since 


f (*i> | z)(z | nw) Mz) = o, 


(5) 


where (z | £l{ip}) is the non-Hamiltonian term in the right-hand side of eq. (4). Evidently, the new 
equation satisfies the homogeneity condition. Moreover, the ’’separability property” also holds: 
the wave function of noninteracting particles is factorized for all times, if it was factorized at the 
initial moment. Nonetheless, the form (4) is not the most general, in [3] we showed that the more 
general form for is 


n{| v>)} = thD 
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where A and r may be arbitrary complex numbers satisfying the restriction A + r = 1. In papers 
[4, 5] we investigated eq. (4) with B = A + , A , B being lowering and rising operators for the 
two-level system and for the harmonic oscillator. It was demonstrated that this equation, written 
in the discrete energy (Fock) basis | n), describes the relaxation to ”pseudothermal stationary 
states, possessing the Planck distribution for the populations of energy levels, but nonzero off- 
diagonal elements of the density matrix. This is not surprising, since eq. (4) relates to pure 
quantum states. However, it is more natural to describe relaxation processes in terms of density 
matrices , in order to investigate the evolution of mixed quantum states. Therefore, our next goal 
is to obtain nonlinear master equations originating from eq. (6) and its special cases. 

Starting from equation (6) and considering the evolution of the pure state density matrix 
(z\ P4,\ z ’) = (z I ip) (ip I z') governed by this equation and replacing | ip)(tp \ by p we obtain 


Wt {z 1 ' 1 z ' ] 


, (£ 1 B*p A* £) 
(z I p I z) (z'\p \ z') 


— r 


VAAiLpA _ <^ K> 

(z'\ p \ z) (z'\p\ z) J 


- D\(z\BAp\z') - D\ m (z\pA+B + \z') - Uz | [fi , p] \ z') , (7) 


which preserves the trace, normalization and hermiticity for an arbitrary initial density matrix. 
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2 Application: The harmonic oscillator in the coordinate 
representation 


This section is devoted to investigating a special case of the general equation derived in Sect. 1. 
We shall limit ourselves to one dimension, identifying the ket- vector | z) =| x) and apply the new 
equation to the harmonic oscillator Hamiltonian 


Li fi 2 d 2 m 2 


- f x 


(8) 


we choose operators A and B as A = B+ = £ + sx, where s,D,r, A, are assumed to be real 
constants. Now, designating the elements of the density matrix in the coordinate representation 
35 p{ x iV) = p*(y,x) and introducing the notation 


H{p } = -~{x\ [H,p] | y) 


% J 

2 —(Pxx(x,y) - p n (x,y)) - ^(V(x) - V(y)) 


( 9 ) 


_ dp _ dp _ d 2 p 
Px = dx ’ Py = dy ’ Pxx = 

we can write eq. (7) as 


( 10 ) 


df l p(x,x) 

, + «y[/>r(y,y) + />»(y,y)] , 

+ 7 r + 2s 

p{y,y) 

, Plxi^y) ± P ly( X ,y)\ . , , / X , 

r P*(x y) J [Parx(^,3/) + Pyy{ x iy)\ ■ (11) 

All derivatives in this equation must be calculated for the independent variables x and y , and only 
after that one should consider x = y in the functions p x ,p y , and p xy . 

For the probability density P(x) = p(x,x) eq. (11) results in the following Fokker-Planck 
equation 

dP dj n (d 2 P d , A 

dt + dx - D { dx 2 + 2,5 dx ^ xP ^\ ’ (12) 

where J(x) = - \p x (x,x) - p y (x,x)\ is the current density, 2 sDx is the drift velocity, and 
D is the diffusion coefficient. Eq. (12) clearly shows that the total probability is conserved in 
time. Eq. (11) admits a special class of solutions in the form of Gaussian exponentials, 


p(x,y) = exp 


-ax + cxy 


^ a'y 2 4- bx + b*y + $ 


(13) 


with c and <I> being real functions of time, while a and b may be, in general, complex. Putting 
expression (13) into eq. (11) we get the equations for the coefficients a(t) and c(<), which do not 
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contain neither the force /(f), nor the function b(t). Therefore, let us assume, for simplicity, that 
f = b = 0. Then, instead of eq. (11) we have the following set of ordinary differential equations 
for the real functions a = Rea, P = Ima, c, and 


a = 2 — ap + 4 D (ac — c 2 — a 2 + s(a — c)} 
m 1 

p = , 2 + - (c 2 - a 2 + /? 2 ) + 4£Z?a/? , 

h m 

c = 2 — f3c — ADac , 

m 

4> - — p + 2 D(c — a + s) 

m 

Here the parameter £ = r — A is introduced . 

The trace of the Gaussian density matrix equals 


rri *• 

1 rp 





consequently, the condition of its conservation is the equation 

1 a — c 


$ - - 


2 a — c 


= 0 


(14) 

(15) 

(16) 
(H) 


(18) 

(19) 


and it is easy to check that this equation is fulfilled. The difficulty of treating the equations for 
the coefficients a, /?, c is connected with their nonlinearity even in the absence of nonlinear terms 
in the master equations (i.e., when D = 0). It is well known, however, that in the latter case the 
equations of motion for average values and the second order moments are linear for any quantum 
system with quadratic Hamiltonian. Therefore we replace the equations for coefficients a, /?,c, by 
the equivalent equations for the variances 

Oxx = (£ 2 ) - <£) 2 1 a vv = ( P 2 ) - ( P ? - = \{xp+ i>x) - {x)(P) 1 ( 2 °) 


these quantities are related to the coefficients of the Gaussian density matrix (13) by 
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(23) 


Another quantity characterizing the quantum state is 


A = 



a pp 



a + c 
a — c 


(24) 
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Its importance is explained by two factors. First, any quantum state must satisfy the generalized 
Robertson - Schrodinger uncertainty relation [6] A > 1. Second, any positively definite density 
operator must satisfy the inequality 


Tr(f?) 

[Tr(f,)Y 


< 1 


(25) 


but for any Gaussian state [7] 


TAP) 

[Tr(p)Y 


A-i/2 


Consequently, the parameter A characterizes the "degree of coherence” of the Gaussian 
A = 1 for pure states, and A > 1 for mixed states. 

The relations inverse to eqs. (21) - (23) read 


(26) 

state: 


so one can check that 
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As a consequence, the equation for the parameter A is 

A - 1 


A = — D(A — 1) 


4s -f 


J XX 
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(28) 

(29) 

(30) 

(31) 


We perceive that A is conserved in time it D = 0. Alternatively, if D ± 0, then A tends to 
the unit value at ( -> oo, provided both coefficients, D and s, are positive. This means that any 
initial state exhibits the relaxation to a pure statel 

The advantage of the equations for the variances is clear: The nonlinear terms are multiplied 
by the diffusion coefficients, which are small in all reasonable situations. Therefore we may use 

the solutions corresponding to D = 0 as the zero approximation, and develop some perturbative 
scheme. 

If Z) — 0, getting rid of a xp and a pp , one arrives at a single third-order equation for the 
coordinate variance, 

d 3 arr. „ da r 


+ 4u 2 


= 0 


whose solution reads 


dt 3 ' dt 

a xx (t) = A + B e 2iwt -f B* e~ 2iwt 


(32) 

(33) 


231 



where A and B are constant coefficient. It is natural to suppose that for D -> 0 the evolution of 
(X xx is given by the same expression, but with the slowly varying in time coefficients. 



<r„(() = A(t) + B+(t) , 

(34) 

where 

B±(t) = B[i) e 2 '"' ± B'(i) e~ 2i “‘ 

(35) 

Now inserting (33) into eq. (28) and neglecting terms of the order of 0(D 2 ) (guessing that the 
n-th derivatives of A(t) and B(t) are proportional to D n ), we can express the covariance a xp as 

follows, 



<7 xp = iu> m B- + 

— (A + Be 2iu>t + B* e~ 2iwt ) + 2 mDs(A + B + ) - mD , 

(36) 


and with the same accuracy we find the momentum variance for eq. (30). 

< 7 pp = (um)\A - B + ) + 8 ium 2 DsB. + 2ium 2 (B e 2i “ l - B* e~ 2iu,t ) 

- ™^[ 2 + <(1+A)] , (37) 

Since in this expression the function A is multiplied by the small parameter D, we must calculate 
it in the zeroth approximation. Thus we get 

A = (^) 2 (A 1 - 4 | B | 2 ) = R 2 , (38) 

being both, R and A, slowly varying functions of time. 

Now, we make the crucial step first proposed in [8, 9]: we put the expression for cr pp in terms 
of A and B into eq. (30) and average both sides with respect to the fast oscillations of frequency 
<jj . Then we arrive at the equation, which contains functions A and R only 

A = 2 To; | -4 <t 0 s + (f - R + “ # Jjp A 

+ (1 + A 2 ) - f) , (39) 

and where <r 0 = h/(2u>m) is the coordinate variance in the oscillator ground state, while T = mD/ti 
is the dimensionless diffusion coefficient. 

Now, averaging eq. (31) we arrive at a closed equation for R(t), 

R = - Tu (R - i) (4^0 s + R - ■ (40) 

Consequently, eq. (39) actually is the first order linear nonuniform equation with respect to 
function A(t), which can be easily solved provided the solution to eq. (40) is known. 

Thus, averaging the equations for the variances over the fast oscillations leads to the effective 
linearization of the equations- governing the evolution of the coefficients of the Gaussian density 
matrix. To get the equation for B(t), one should first multiply both sides of the equation for a pp , 
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by exp(— 2iut), and then average over the fast oscillations. We confine ourselves, however, to the 
evolution of functions A (t) and A(t) (the absolute value of B(t ) can be extracted from eq. (38)). 
It is worth to mention that A(t) is nothing but the energy of quantum fluctuations (up to the 
constant factor and small corrections of the order of 0(D )): 

SE{t) = + 2~ <7 pp = • (41) 

Now if we consider a pure state, R — 1, then the solution for eq. (39) is 

,4(f) = <7 0 + [A(0) - (7 0 ] exp [- ATu>(2a 0 s - f )f] . (42) 

We see that the energy of fluctuations can increase or decrease, depending on the sign of the term 
in the exponential: 

1) If 2oo s — £ > 0, we verify that ,4(oo) = <7o , 6E(oo) = %u/2 and B(oo) = 0, therefore 
asymptoticaly a xx and <j pp do not oscillate. 

2) If 2<7 o'S — < 0 and since reasonably f\j <C 1, A(/) becomes a slowly increasing function with 

time (the same for 6E(t)). From eq. (38), asymtoticaly we have A{t) ~ 2 B(t) <Tq. If B(t) is real 
at t — 0 then it is real for any t > 0, so from eqs. (34) and (35) cr xx (t) = A(t) + 2B(t) cos(2ut) 
or A(t) — 2 B(t) < o xx < A{t) + 2 B(t). But from eq. (38) A(t) — 2B(t) = o-Q/(A(t) + 2 B(t)) ~ 
<rl/(2A(t)). Therefore asymtoticaly a rx oscillates between two values, 

2 

'2X[t) ” <Jxx " ’ (^) 

thence, as time goes on the solution of the nonlinear Schrodinger equation becomes highly squee- 
zed. Note that a similar behavior of the variances (i.e., an exponential increase of the squeezing 
coefficient) is observed in the case of the usual Schrodinger equation for the parametrically excited 
oscillator, when its frequency changes in time [7]. However, in the present case all the coefficients 
in the generalized Schrodinger equation with the functional (6) (or its master equation counterpart 
eq.(7)) do not depend on time, and the increase of fluctuations is caused by the nonlinear terms. 

3 Amplifying — Dissipative Hamiltonians 

Here we shall consider the theory developed many years ago by P. Havas [10], which is quite 
suited to construct Hamiltonians that take into account dissipation and apply it to two examples, 
confining ourselves to the one-dimensional case. 

3.1 The Bateman-Caldirola-Kanai (BCK) Hamiltonian 

The harmonic oscilator with an exponential time- dependent mass m(t ) = m 0 exp(^), is known 
as the BCK Hamiltonian [11] and we shall introduce the phenomenon of friction in it. According 
to [10] when friction is present the Lagrangian that describes the motion, ( that reproduces the 
classical equation of motion) is 

L{q,<r,i) = (1™(0? 2 - + <?^(0) exp[7(f)] . (44) 
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Here V(q;t ) = |m(<)u)^ q 2 and F(t) is a time-dependent external force. The exponential factor 
introduces the friction effects (also external) on the system and for the specific Lagrangian one 
has 

7(0 = y (1 ~ exp(-<0) . ( 45 ) 

as lim t _ >00 7(<) = 70/C and lim c _o7 (f) = 70 1, so, recovering previous results [12]. We shall 
consider F(t ) as being a stochastic force whose mean in an ensemble is null, ( F(t )) = 0, and that 
its correlation is Markovian , ( F(t)F(t ')) = 2 d8{t — t'). 

According to the quantization procedure of [12] one obtains the system Hamiltonian 

A - 

^ exp{-(t)+ \m 0 u> 2 Q 2 exp((t) +^{Q,P}exp(-(t) + QF(t)exp(7{t)/2) (46) 
Zmo Z 4 

where P,Q are canonicaly conjugated operators, [Q,P] = ih. The physical position and momen- 
tum operators are related to those through 

q P h V s = $exp(-7(*)/2) , Pphy, = •f >ex p(— 7(0/2) ( 47 ) 

In order to obtain the equations of motion for the operators P and Q in the Heisenberg pic- 
ture we first do an unitary transformation in Schrodinger equation with operator S(Ct/2) = 


exp(— P})i which leads to the new Hamiltonian 

Tt 2 \ 1 

/C(M;0 = 7 — + X m oWo ^ 2 + 7 ( 7 oexp(-C 0 + C){^,^} 

Zttiq Z 4 

+ xF(t)exp ^(7(0 - CO] ( 48 ) 

with 

q P h v , = ^exp [- ^(7(0 + CO] 

Pphy, = ffexp[-^(7(<)-C0] • ( 49 ) 

The equations of motion for x(t ) and H(t) are solved assuming xjj = u{t)x 0 + v(t)H 0 + w (t)> 
leading to the set of linear differential equations 

^7 - H* + yexp(-2(t) u = 0 (50) 

^7 - + Y ex P(“ 2 C0 v = 0 (51) 

^ - Ho + ^exp(-2C0 w = “ ^^ ex P f^(7(0 “ CO] ( 52 ) 

dt 2 l J m 0 J 

where = £ — u>q and with the following initial conditions 

u(0) = 1 and u(0) = -(70 + C) (53) 

w(0) - 0 and u(0) = — (54) 

TTIq 

w(0) = u>(0) = 0 . (55) 
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The exact solutions of the above differential equations are : 

” w = [(" _ i" 5 ) K -(i) + exp (- c ‘ 

“ [(" " 2 ? - jM!) - i 1 -' (|)] K -(f exp(-co) 


,(f) = - K -(!M! exp( - c( 0 ] 

’ (t) = nw! fo exp {j['>' (( > )_ ^ 1 ]} [ I ‘'(^ ex P(-CO) K >-(Sexp(-C«i)) 


- K 


- (2^ ex P(-C^i))] F(ti) dti 


where Ij,(.) and K 1/ (.) are the modified Bessel functions of the first and third kind, respectively 
[13]. The parameter u is 

^ if 2 > 

" = ' 0 if § = w 0 . (59) 

, if ^ < ujo 

Similarly the momentum operator in the Heisenberg picture is Ttn{t) — /x(t)x 0 + j3(t)j r 0 + £(t), 
where 


H(t) = m 0 ii(t ) - ^ [70 exp(-C<) + (] u(t) 

0 (t) = m 0 v(t) - ^ [ 7o exp(-Ct) + C] v(t) 

f (0 = m 0 MO - ~ [70 exp(-<<) + C]w(< 


Considering that the energy increasing prevails over energy dissipation, 7o /C <C 1, one has the 
approximative solutions, 


xh {0 — ^cosh(fl 0 0 + ^°20 ~ sinh(fl 0 t)^xo + (^ — ^ sinh(fl o 0^o 

“ exp ( 1 2 ^) l exp (“ 2 ^) ( 61 ) 

*»(*) = [-^ £ sinh(Qo<) + yp(cosh (n 0 t) - ^-sinh(fV))] x 0 

+ (cosh(fV) - sinh(H 0 t))iro 

+ exp(^l) £ ex P^~y~ ^ sinh[fto(«-<i)] - cosh[f) 0 (< — <i)]J F{t\)dt\ (62) 
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(63) 


and the variances for q p h ya and p p h y , become 


Um(A Xl#) 

~ — - — expf- 

160 0 

-K- 2 n.M{a(^ ; - i )«i h ( 

hu 0 \ i 
2 k B TJ 


exp 

<-?) 

( 2 ^ + 1 )(' + j) 

l }-0 

0 1 

>0 

lim(A Xl„) 

t — ►CO 

c 

160, 

-exp[(C + 2n„)(|{^(^- 

+ 1 ) COth ( 2 (c b t) + 


exp( 


) 2 n„ ‘)( 1 + t)j 


' c 

v 2 flo 

> 0 

verifying squeezing even 

in occurence of 

a weak dissipation. 





(64) 


3.2 The optical parametric oscillator (OPO) 

Using the same method we can treat the OPO introducing in the Hamiltonian the dissipation of 
the cavity, hence, 

H{A + , A, t ) := fi(t)A + A + {f 2 {t)(A + ) 2 + f 3 {t)A+ + h.c .} (65) 

where the mathematical operators A+ and A are related to the physical creation and destruction 
operators by a = A exp (-At/2) and a + = A + exp (-Xt/2). Moreover, /i(t) = lo 0 is the 
mode frequency in the cavity, = Kexp(— 2 iuit) + iA/4, where k and u> are, respectively, 

the intensity and the frequency of the pumping field, while A is the damping constant of the 
cavity. f 3 (t) = F(t) e Xi /(2w 0 ) 1/2 • The force F(t) is assumed to be a Markovian stochastic force: 
< F(t) > = 0, < F(t)F m (t') > = < F(t)*F(t') > = 2 dS(t - t’) and < F(t)F(t') > - 0, and the 
parameter d is related to the temperature of the cavity, 



Considering the pumping at resonance, u>o = u?, and for k\/u)q 1, the solution of Heisenberg 
equations for operator An{t) is 

A H {t) = u(t)A(t) + v(t)A + (t) + w(t) (67) 


where 

u(t) 

v{t) 

w{t) 


e ,u/ 0 < cosh( 2 nt) + |e 
t sinh( 2 /c<) (cosh 2 7 e 

— cosh (7 — 2 Kt) 

Jo 


~ twot (sinh 2 7 cosh( 2 /c t) — 1 sinh 7 Cosh 7 sinh( 2 «t)) — 
-»w 0 t _ s i n fi 2 e ,w °^ — i sinh( 27 ) sinh( 2 «<) cos( 2 a;ot) 

f 3 (t')exp[-iu 0 (t - t')]dt' - c.c.l 


c.c. 


( 68 ) 

(69) 

(70) 


and tanh( 27 ) = A/(2u>o). 
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With the above solution we present the asymptotic ( t — ► oo) mean values of several quantities, 
after taking the average over the high frequency oscillations of the field (cj 0 ) and considering 
4k/ A <1: 

1) Energy, 

hm (E) = u 0 e- Xt (A^A H U ~ j —(l K /\) 2 ( 71 ) 

2) Variances of the two quadratures, 


\im{AX h2 ) t ~ 


_d 1 

2Xu>o 1 — (4k/A) 2 


(72) 


Although the asymptotic values of the variances of the quadratures are both above the value 
1/2, at initial times squeeze is seen in the fast oscillations (before averaging over the mode fre- 
quency). In eqs. (71) and (72) one verifies the effects of the dissipation-amplification process 
through the quantity 4k/A: 

i) If 4k/ A < 1 (very weak pumping compared with the cavity dissipation), the thermalization, 
represented by the parameter d , dominates in the physical expressions at equilibrium. 

ii) On the other side, when 4 k/A is close to 1 (strong pumping) the factor (1 - (4 k/A) 2 ) -1 dom- 
inates the strejights of the asymptotic values, increasing dramatically the energy and fluctuations, 
as is expected to occur at resonance. 


4 Summary 

We have presented two different formulations of quantum dynamical equations that show squeezing 
in the variances of the conjugate canonical operators. In the first one we considered a general- 
ization of the Doebner-Goldin nonlinear extension of the Schrodinger equation and we verified 
that although the parameters that enter the nonlinear part of the equation are constant in time, 
squeezing occurs, essentially due to the nonlinearity. Moreover, the master equation shows the 
surprising feature that any initial mixed state relax to a pure state! 

In the other approach we introduced the dissipation phenomenon into the Hamiltonian for- 
malism by starting with a conveniently defined Lagrangian, as proposed by P. Havas [10]. We 
considered two familiar time-dependent Hamiltonians, the BCK and the OPO. The BCK Hamilto- 
nian has a time-dependent mass and it displays amplification of energy and squeezing of variance 
of momentum or of position, although uncertainty is preserved. The dissipation was introduced 
and the effects are seen in eqs. (63)-(64). 

The second Hamiltonian is the OPO, describing a single mode in an electromagnetic cavity 
with pumping at resonance. Dissipation is introduced to take into account the loss in the cavity 
walls. As an expressive result we verify that the asymtotic physical expressions depend, essentially, 
on the factor 4 k/A, k representing the pumping and A, the dissipation. 
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